Low Energy Chiral Lagrangian in Curved Space-Time from the Spectral Quark Model 
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We analyze the recently proposed Spectral Quark Model in the light of Chiral Perturbation Theory 
in curved space-time. In particular, we calculate the chiral coefficients Li, . . . ,I/io, as well as the 
coefficients Ln, L12, and L13, appearing when the model is coupled to gravity. The analysis is 
carried for the SU(3) case. We analyze the pattern of chiral symmetry breaking as well as elaborate 
on the fulfillment of anomalies. Matching the model results to resonance meson exchange yields 
the relation between the masses of the scalar, tensor and vector mesons, 

4-y/3/iVc7r/,r. Finally, the large- A^c limit suggests the dual relations in the vector and scalar channels, 
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I. INTRODUCTION 

The low-energy structure of QCD in the presence of 
external electroweak and gravitational sources is best de- 
scribed by Chiral Perturbation Theory (ChPT) X,^^ 
0,0 (for review see, e.g., Ref. @). In the meson sector, 
the spontaneous breaking of chiral symmetry dominates 
at low energies and systematic calculations of the corre- 
sponding low-energy constants (LEC's) have been carried 
out in the recent past up to two loop accurac y [3,H,|^ ll0 | . 
or by using the Roy equations (see also |l2lll3j |). For 
strong and electroweak processes involving pseudoscalar 
mesons the bulk of the LEC's is saturated in terms of 
resonance exchanges jl4j . which can be justified in the 
large- A'^c limit in a certain low-energy approximation [Tsj l 
by imposing the QCD short-distance constraints. In the 
case of gravitational processes similar ideas apply |^ , al- 
though less information is known Ig*]. Nowadays, ChPT 
can be used as a qualitative and quantitative test to any 
model of low-energy hadron structure. 

In the quest to understand the microscopic dynamics 
underlying the LEC's, their calculation in chiral quark 
models has been undertaken many times 0, 0, 0, [23, 
El 111 111 111 111 111 mm m 113. The effort has 
been made to compute Li, . . .Lio, which correspond to 
the flat-space-time case. The calculation of Ln, L12 and 
L13, encoding the coupHng to gravitational sources, has 
not been considered so far. Roughly speaking, these cal- 
culations are generally described in terms of some long- 
wavelength expansion of the fermion determinant asso- 
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ciated to the constituent-quark degrees of freedom. A 
detailed scrutiny shows, however, that the implementa- 
tion of the necessary regularization is not always satis- 
factory from several viewpoints. The regularization of a 
low-energy chiral quark model corresponds to a physical 
suppression of the high-energy quark states. This can 
be achieved in a number of different ways, e.g. by cut- 
offs, form factors, or momentum-dependent masses, pro- 
vided they do not break symmetries such as the gauge 
invariance and chiral symmetry. Thus, the regulariza- 
tion should not be removed in the end. In such a situa- 
tion, where the high-energy quark states are suppressed 
above a certain scale A, one should expect a power-like 
behaviour A"/(5" for any large-momentum external leg 
of the quark loop in the high-momentum limit. In the 
language of the parton model this high-energy behaviour 
corresponds to the onset of scaling. 

As a matter of fact, one of the questions which could 
not be answered by low-energy calculations concerns the 
low-energy resolution scale where these models are sup- 
posedly defined. Actually, in order to properly answer 
this question one should look instead into high-energy 
processes and demand parton-model relations on the con- 
stituent quarks. As pointed out in Ref. 31], a sensible 
scheme is obtained by demanding that the momentum 
fraction carried by the valence quarks in a hadron sat- 
urates the energy-momentum sum rule. Once this ini- 
tial scale is defined one can use the QCD evolution to 
compute an observable at a higher scale. This way the 
QCD radiative corrections are incorporated. In fact, us- 
ing the analysis of the Durham group carried out a decade 
ago (32] for the case of the pion, one obtains the result 
that the valence quarks saturate the energy-momentum 
sum rule at fio = 313 MeV if the LO DGLAP QCD 
perturbative analysis is carried out. Although this scale 
looks quite low, the impressive agreement obtained for 
the parton distribution functions of the pion after the 
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DGLAP evolution in LO [33, |3J and NLO |3J (see also 
Ref. i3LL and Ref. [S^ where the comparison to the E615 
data [2a| is made) supports this interpretation of the low 
resolution scale. Moreover, using that scale, the pion 
distribution amplitude j37il and the off-forward general- 
ized parton functions [23 agree well also with the recent 
transverse lattice calculations [s^ [s^ , which presumably 
incorporate the non-perturbative evolution. 

A proper identification of the low-energy matrix ele- 
ments entering the high-energy processes is grounded on 
the absence of logarithmic corrections in the low-energy 
model in the high-energy limit, since the proper QCD 
radiative logarithmic corrections are automatically and 
completely incorporated by the QCD evolution. Not sur- 
prisingly, this condition imposes severe constraints on the 
kind of admissible regularization schemes. In a recent 
work the Spectral Quark Model (SQM) has been pro- 
posed [i^, llj , implementing the so-called spectral regu- 
larization (see below) complying to these power-like short 
distance constraints. 

In the present paper we extend the SQM to the SU(3) 
flavor group and include finite current quark mass. In- 
stead of using the construction of vertices based on the 
Ward-Takahashi identities, employed in Ref. [iol |4H. it 
is by far more convenient to define the effective action 
depending on the non-linear pseudoscalar meson fields 
in the presence of external scalar, pseudoscalar, vec- 
tor, axial, and gravitational sources. The latter have 
never been considered in chiral quark model calcula- 
tions. This effective action is defined in Sect. ^ We 
also show in Sect. IIIII how one can explicitly eliminate 
the spectral function in terms of the quark momentum- 
dependent mass and wave-function renormalization. Fol- 
lowing the standard procedure we perform the gradient 
expansion of the spectral-regularized fermion determi- 
nant both for the anomalous (Sect. IIV|I as well as the 
non-anomalous sectors in curved space-time (Sect. Iv|l . 
As a consequence the structure of the energy-momentum 
tensor may be analyzed. Remarkably, our spectral regu- 
larization method complies to the QCD anomaly with- 
out removing the regularization. Therefore the stan- 
dard Wess-Zumino-Witten ,4^ ,4^] term is generated for 
a finite regularization. In the non-anomalous sector we 
find, through the comparison to the standard chiral La- 
grangian [30], that the low-energy constants at 0{p'^) 
associated to terms which are non-vanishing in the chi- 
ral limit are completely independent of the regularization 
details. The LEC's associated to terms carrying the cur- 
rent quark mass coefficients do depend on the particular 
ansatz for the spectral regularization, and we evaluate 
them using the regularization based on the meson dom- 
inance of form factors 41]. Such a model has provided 
a satisfactory description of the quark self-energy of the 
recent lattice data [ij]. Finally (Sect. IVIljl . we also con- 
front the large-TVc relations [l^ and discuss the conse- 
quences of extending the present model to include these 
constraints. The Appendix IXI contains details of the for- 
malism in the curved space-time. 



II. THE EFFECTIVE ACTION OF THE 
SPECTRAL QUARK MODEL 

In a recent work the spectral quark model (SQM) has 
been introduced |40ll4l| . The approach is similar in spirit 
to the model of Efimov and Ivanov jl^ , proposed many 
years ago. It is based on the formal introduction of the 
generalized Lehmann representation for the quark prop- 
agator. 
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where p(ui) is a (generally complex) quark spectral func- 
tion and C denotes a suitable contour in the complex lo 
plane. The function M(p'^) is the quark self-energy, while 
Zip^) is the quark wave-function renormalization. In the 
case of analytic confinement, i.e., when the propagator 
does not have poles, a sensible definition of a constituent 
quark mass is (from now on wc drop the index C from the 
m integral, which is implicitly understood to run along 
the contour C) 
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As discussed at length in Ref. Al} , the proper normaliza- 
tion and the conditions of finiteness of hadronic observ- 
ables are achieved by requesting an infinite set of spectral 
conditions for the moments of the quark spectral function 
p{lo), namely 

Po = J dLjp{uj) = 1, (3) 

Pn = J dcjuj^piui) = 0, for n = 1,2,3,... (4) 

Physical observables are proportional the zeroth and the 
inverse moments, 

p-k = J dLJUj-''p{uj), for fc = 0,1,2,3,..., (5) 

as well as to the "log moments", 

p'^^ I d^log(c.V^2)^n^(^) 

duj\og{uo^)uo''p{uo), for n = 1,2,3,4, ... (6) 



Obviously, when an observable is proportional to the di- 
mensionless zeroth moment, po = \, the result does not 
depend explicitly on the regularization. The spectral con- 
ditions ^ remove the dependence on the scale in 10, 
thus guaranteeing the absence of any dimensional trans- 
mutation. The only exception is the Ot/i-log moment. 



Po(m ) = / dwlog(w /p, )piuj). 



(7) 
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which does depend on a scale fj, and is not regularized by 
the spectral method (see the discussion below) . No stan- 
dard requirement of positivity for the spectral strength, 
p{lo), is made. Unlike other regularizations, such as the 
dimensional regularization or the ^-function regulariza- 
tion, the spectral regularization is physical in the sense 
that it provides a high-energy suppression in one-quark- 
loop amplitudes and is not removed at the end of the 
calculation. It also improves on a Pauli-Villars regular- 
ization, because it complies to the factorization property 
of correlation functions, form factors, eic, in the high- 
energy limit, i.e., it guarantees the absence of logarithmic 
corrections to form factors. The phenomenological suc- 
cess of the SQM in describing structure functions of the 
pion, generalized parton distributions '35', '46], and the 
pion light-cone wave function 37i j4^] suggests that the 
whole scheme deserves to be thoroughly pursued further. 

In Ref. it was argued that there are a number 
of terms in the one-quark-loop effective low-energy chi- 
ral Lagrangean which correspond to taking the infinite- 
cut-off limit. The terms with explicit chiral-symmetry 
breaking do not correspond to this class. The purpose of 
this paper is to analyze these terms, which are specific 
both to the regularization and the choice of couplings 
in the spectral quark Lagrangian. For completeness we 
also consider the gauge couplings and gravitational cou- 
plings, which allows us a determination of all low-energy 
constants in the SU(3) sector in the SQM approach. 

The effective action complying to the solution of the 
Ward-Takahashi identities via the gauge technique of 
Delbourgo and West 0] corresponds in our case to 
the minimum substitution prescription for the spectral 
quark. It yields a quark fermionic determinant of the 
form 



T[U, s,p,v,a, g] — —iNc / diLjp{uj)Ti- log (iD) 



(8) 



where the Dirac operator is given by 

iJD = — LoU^ — rho + (J) + (^75 — s — i^^p) 



(9) 



The derivative is frame (local Lorentz) and general- 
coordinate covariant and it includes the spin connection 
(see Appendix 1X1 for notation). The symbols s, p, w^, 
and denote the external scalar, pseudoscalar, vector, 
and axial flavour sources, respectively, given in terms of 
the generator of the flavour SU(3) group, 



a=0 



Aa 

' 2 ' 



(10) 



with Aa representing the Gell-Mann matrices. The tensor 
gfj^u is the metric external source representing the cou- 
pling to a gravitational field. The matrix = If^, and 



U 



tV2'S>/f flavour matrix representing the 



pseudoscalar octet of mesons in the non-linear represen- 
tation. 
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(11) 



The matrix ttiq = Diag(m„, m^, ms) is the current quark 
mass matrix and / denotes the pion weak-decay constant 
in the chiral limit, to be determined later on from the 
proper normalization condition of the pseudoscalar fields. 
For a bilocal (Dirac and flavour matrix valued) operator 
A{x, x') one has 



TrA 



-gii{A{x,x)) , 



(12) 



with tr denoting the Dirac trace and ( ) the flavour trace. 
Moreover, g = detg^^ is the determinant of the curved 
space-time metric. Finally, in the second line of Eq. 
we have introduced the Dirac operator D correspond- 
ing to the external fields only. The Ua{^) is taken into 
account by extending the matrix to the U(3) sector, 
U ^ U = [/e'''8/(3/) with detU = 1, adding the cus- 
tomary term 



C 



P 2 



61- ^ [log detC/- log detC/t] I . (13) 



The Dirac operator given by Eq. @ transforms co- 
variantly under local chiral transformations (see Ap- 
pendix . 

Formally, in the flat space-time the effective action ((H)) 
looks quite familiar and we should point out here that 
the main difference with similar actions, such as, e.g., 
the one of Ref. |23| , is related to the reg ularization proce- 
dure. Actually, the method of Ref. 21] consists of taking 
p{uj) = 6{lo — Mq) with Mq being the constituent quark 
mass. This choice satisfies the normalization condition 
Pa — 1, but does not comply to the Pn = spectral 
requirements. The problem can be avoided if one uses 
suitable regularization methods, such as the dimensional 
or C-function regularization, but then logarithmic correc- 
tions to form factors are generated and the well-known 
Landau instability found long ago in Refs. Isoj ) sets 



The pion form factor obtained from the ^-function reg- 
ularization used for instance in Ref. 21] for t = —Q^ 
becomes, in the chiral limit. 



X{1 - x)Q2 + M,2 



Q 



p^ 



where the pion weak-decay constant is given by = 
4:NcM^ log(/i2/M|)/(47r)2. While the proper normafiza- 
tion -F(O) = 1 is obtained, at large momenta one has a 
logarithmic behavior, F{Q^) log((5^), instead of the 
power-like behaviour, which poses a problem. On the 
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other hand, the spectral regularization method yields |43 

logs present. This twist expansion property allows us to 
extract in a clean way the low-e nerg y matrix elements 
relevant for high-energy processes [31|. 

Given the fact that the integration contour is in general 
complex, passing to the Euclidean space and separating 
the action into the real and imaginary parts becomes a bit 
inconvenient. Instead, we take the full advantage of the 
Minkowski space and introduce the auxiliary operator, 

- iD5 = 75 - ^U^^ - Too + ^ - 75^ - s + 175^) 75, 

(14) 

which corresponds to the Hermitean conjugation in the 
Euclidean space. Thus, the normal parity action is given 

by 



Nc / dwp(w)Trlog(DD5 



(15) 



III. RELATION OF SPECTRAL MOMENTS TO 
QUARK MASS AND NORMALIZATION 

A potential disadvantage of the spectral regularization 
is that the inverse problem, i.e. the problem of finding 
the spectral function p{uj) from the known moments, does 
not always have an easy explicit solution or perhaps has 
no solution at all. In this section we show how the nega- 
tive moments and the log-moments can be translated into 
the integrals involving the quark mass function, Af (p^), 
and the quark wave function renormalization, Z{p^). Let 
us start with Eq. ^ and assume that the set of spectral 
conditions is met 



dww'V(tj) = (5„o, n = 0, 1,. 



(16) 



Then, the following identity, proved by induction, holds 



dio- 



Rationalizing the denominators yields 



(17) 



— U!^ 



(18) 



We have two cases of odd and even n. For n — 2k we 
find 



9 ? 



p2 - M{p'^Y 



Defining 



? 2 
— UJ^ 



i>p^^-\ 

(19) 
(20) 



and comparing coefficients of powers of p in Eq. refrela 
produces the identities 



L2k{p^) = P^'Zip^) 
L2k+i{p^) - P^'Zip') 



1 



„2fc-2 



p2 - M{p^Y 



p2 - M{p- 



,2\2 ■ 



(21) 



The case n — 2k + 1 produces the same relations. 

The following recursion relations follow directly from 
the spectral conditions 



w"p(w) 

p2 _ ^2 



p'^ — 



n > 2, 



(22) 



which are obvious when on the right-hand side we write 
j)2 — {jp' — up) ^ up' . We now pass to the Euclidean space, 
^2 _ p2 _^ — and get 

/•OO 

da>w"log(c^2)^(^) ^ j dplLnhpl) (23) 



/ dp\p%Ln-2{~p]s)- 

Jo 



Thus, we have obtained the log-moments in terms of Z 
and M . The negative moments are simply derivatives of 
the quark propagator at the origin, 



n-l 



p=0 



n = l,2,...(24) 



The derivative is computed taking p^ — fj). Thus, given 
the quark propagator S(j)) we may just use formulas 
(|23I24|I to translate negative moments and log-moments 
without ever having to specify explicitly the spectral 
function. This is a rather remarkable feature of the spec- 
tral approach. The expressions for (gg) (the quark 
condensate for a single flavor), and B (the vacuum en- 
ergy density) in the chiral limit are 



f = 

{m) = 

-B = 



AN, 
(47r)2 

(47r)2 
74^ 



duiuP p{ijj){'- logcj^), 
dujuj^ p{uj){— logcj^), 
dujuj'^ p{uj){— logcj^) 



(25) 
(26) 
W,(27) 



respectively. Here, 9^^ is the energy momentum tensor 
(see also Sect. IV We get, for instance. 



dpi 



or 



ANc 

jAnfj, 



m{-pIY^pI{z{-pI)-i) 

pl + M{-pir 
Zi-Pl)M{~pl) 



(28) 



(29) 
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In Eq. (|29(l we recognize the usual formula for the quark 
condensate found in non-local models. On the other 
hand, Eq. (|28ll is different from analogous quark-model 
expressions |5ll|52j. The reason is that, strictly speak- 
ing, the above formulas should only be used for functions 
M(p^) and Z{p^) complying to the generalized Lehmann 
representation, Eq. Q , with the spectral density satisfy- 
ing the spectral conditions. 

One can use similar manipulations to get the pion elec- 
tromagnetic form factor obtained in Ref. l4l|| . For space- 
like momentum, ~ —Q^, we obtain 



with 



where 



dx 



(30) 



■ ^^2^M{-Plf-Pl{Z{-Pl)-l) 



Pl + M{-PlY 



p2 



(31) 



Note that the inversion procedure used in Ref. (41| to 
determine the spectral density from vector meson dom- 
inance (the Meson Dominance version of the SQM) is 
linear, whereas written in terms of M and Z becomes 
highly non-linear. 



IV. CHIRAL ANOMALIES 

One of the major advantages of the spectral regular- 
ization is that it makes hadronic observables finite and 
scale independent, a basic requirement of any regular- 
ization procedure. However, that does not necessarily 
mean nor imply that the full effective action in the pres- 
ence of external fields is finite, since even in the case of 
the vanishing pion fields, = 1, we have non- hadronic 
processes. Actually, it turns out that the photon wave 
function renormalization [Sj is proportional to pg, thus 
it depends on the scale fi and therefore diverges in some 
regularization schemes (such as the dimensional regular- 
ization). This scale dependence arises also in other non- 
hadronic terms of the effective action. 

In Ref. it was checked that the 7r° 27 and 7 
Stt decays comply to the correct values expected from 
the chiral QCD anomaly. With the help of the effective 
action, Eq. ||SJ), we now want to show that this is also 
true for all anomalous processes. In order to understand 
the role of regularization, it is instructive to compute 
the chiral anomaly first. Next, we will show that in the 
presence of external fields the anomaly does not depend 
on the pion field U and thus coincides with the anomaly 
in QCD due to the spectral conditions pi = P2 = P3 = 

P4 = 0. 

Under chiral (vector and axial) local transformations 
the Dirac operator transforms as 



ev{x) ^Y.''v{x')^a, eA(x) =^el(a;)Aa. (33) 



Infinitesimally, we have 



(34) 



If we make a chiral transformation of the effective ac- 
tion © without any additional regularization, we get 



5S = -iNc^r j dLup{uj) [SDT)-^] 



(35) 



If we assume the cyclic property of the functional trace 
we get a contribution from the axial variation only, 

6aS = Aa = J d^xti J duj p{u)) {2iaj5) 

= po / d'^xti-{2iaj5), (36) 



a result which, due to the infinite dimensional trace |53 . 
Is^ . is ambiguous even in the presence of the spectral 
regularization. Thus, to get rid of the ambiguity we have 
to introduce an extra regularization. As is well known, 
there is no regularization preserving the chiral symmetry, 
thus the anomaly is generated. 

The calculation can be done by standard methods. A 
very convenient one is the C- function regularization 5^, 
which computes the anomaly directly in terms of the 
Dirac operator itself (and not its square) and does not re- 
quire any redefinition of the Dirac 75 matrix. This yields 
the equation 

SaS = Aa = Tt J dcjpiuj) (2^75 [iD]°) (37) 
= / d'^xiti / dujp{uj){2ia{x)j5{x\r>"\x)), 



where the zeroth power of the Dirac operator is under- 
stood as an analytical continuation which can be written 
in terms of the Seeley-DeWitt coefficients for the Dirac 
operators [ssj . 



1D"| 



D 



3+i£v(a;)-ieA(a;)76£)g-iev(a;)-ieA(a;)75 



(32) 



^ (3^{^^^' + ^(D'r^ + r-°'r, + r^D^) 

+ 1 + {r^T^f + T^TlT,) }, (38) 

where = -^{7^, D} and the operator D acts to the left. 
The result for general couplings in four dimensions has 
been obtained from Ref. [53 ■ Direct inspection shows 
that since the (jj-dependence is given by iT) = iD — ojU^, 
the result can be written as a sum of an w-independent 
term and a polynomial remainder, 

Aa = J dujp{u}){AA[v,a,s,p]+AA[v,a,s,p,uj,U]) 
= poAA[v,a,s,p], (39) 
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where the w-dependent polynomial term vanishes due to 
the spectral conditions. This shows that the anomaly 
of the spectral quark mode coincides with the anomaly 
of QCD after introducing an additional suitable regular- 
ization, regardless of the details of the spectral function. 
This result is common also to nonlocal models when one 
evaluates anomalies [s^ l57| . This is an important point 
since if the effective action r[U, s,p,v,a] in Eq. ||HJ) is both 
chiral symmetric and finite there is apparently no reason 
for anomalies. We will see below how and where these 
divergences arise. 

To see now how the standard Wess-Zumino-Witten 
11^ term arises in the present context, let us consider for 
simplicity the chiral limit toq = and set the external 
fields to zero and work in flat space, so that iD = i^. A 
convenient representation can be obtained by introducing 
the field 



5 _ „JtV275*// 



(40) 



Higher order corrections to the abnormal parity com- 
ponent of the action involve negative spectral moments. 
For instance, the terms 0{p^) and higher are regularized, 
and involve p-2 for terms with no quark mass terms and 
for terms containing one quark mass. This is in con- 
trast to the approach of Ref. where the infinite cut-off 
limit is considered for a constant constituent quark mass. 
In this regard let us also note that for the unregularized 
abnormal parity action one would get the transition form 
factor 



8AB 



(47r)2/^ 



dx 



1 



(1 - x)xQ^ + 2Af 2 



which satisfies the proper anomaly condition, ^",^77* (0) = 
1 / {iiT^ f.^) . Again, a log dependent term is obtained 
at high virtualities (see also Ref. [sJl), in contrast to 
the correct twist expansion generated by the spectral 
method [iHl. 



interpolating between the vacuum 
matrix C/p_ 



uL 



1 , and the full 



— U . Then, we have the trivial but useful 
identity for the vacuum-subtracted action. 



LOW ENERGY CHIRAL EXPANSION OF 
THE ACTION 



r[u,s,. 



..] - r[i,s,...]= (41) 

- i Na [ dt^ [ dcj p{lj)Tt log {iD - LuU^^) 
Jo Jc 



' dt / da;/9(aj)Tr 
Jc 



dU^ 



1 



dt W-ujU^ 



Using the representation in Eq. (|42|) and the formulas of 
Appendix^the result can be obtained straightforwardly. 
Since we are interested in abnormal parity processes it is 
enough to identify the terms containing the Levi-Civita 
tensor e^^ap, which due to the Lorentz invariance re- 
quires at least four derivatives. Taking into account the 
fact that the derivative operator acts to the right we get 



-(4) 



-iNr I dt dhjpiuj) 
'0 Jc 



d^x 



d^k 



(27r)4 [fc2 



Tr <j -Lo-i^ul^LO \LoUli^Ui^ 



,215 



(42) 



After computation of the traces and integrals we finally 
find 



-(4) 
ab 



Po- 



dt 



,4 



(43) 



X {ul'^uld^^UtUld^UtUld'^UtUld^Ut), 



which coincides with the WZW term if the spectral nor- 
malization condition pQ — 1 is used. External fields can 
be included again through the use of Eq. yielding 
the gauged WZW term in the Bardeen subtracted form. 
Actually, the difference r[[/, s,p, a] — r[l,s,p, z;,a] is 
finite and preserves gauge invariance but breaks chiral 
symmetry generating the anomaly of Eq. H39|l . 



The chiral expansion of the action, Eq. © , correponds 
to a counting where the pseudoscalar field U and the 
curved space-time metric g^'^ are zeroth order, the vec- 
tor and axial fields and are first order, and any 
derivative 9^ first order. The external scalar and pseu- 
doscalar fields s and p and the current mass matrix mo 
are taken to be second order. In chiral quark models 
at the one loop level this chiral expansion corresponds 
to a derivative expansion. With the help of the action 
of Eq. |(SJ) one can compute the derivative expansion in 
curved space-time (see Appendix IXI for details). 



S ■ 



-gC{x) 



(44) 



where the effective chiral Lagrangian in the Gasser, 
Leutwyler and Donoghue form |j, |^ reads 



£ = /:(0)+£(2.9)+£(2,fl)+£(4,<,)^^(4,fl) 



(45) 



with the metric (upperscript g) and curvature (upper- 
script R) terms explicitly separated. The zeroth order 
vacuum contribution reads 



^ - (47r)2''4' 

where the vacuum constant is given by Eq. H27|) . 



A. Metric contributions 

The metric contributions read 



(46) 



£(2,9) 



(D^U^D^U+ix^U + U^)), (47) 
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and 

£(4,9) ^ L.iD^U^D'^U)^ + L2{D^U^D,Uf 
+ L:i{{D^U^D,uf) 
+ L4{D^U^D^U){x^U + U^x) 
+ L5{D^U^D^U{x^U + U^x)) 
+ U{x^U + U^x? 

+ Lr{x^U~U^xf+Ls{{x^Uf + {U^xf) 
- iLg{F/;:^Df,UDM^ + F^.D^U^DM) 

+ miiF/^^f + (F^^f) + H,{xh)- (48) 

We have introduced the standard chiral covariant deriva- 
tives and gauge field strength tensors, 



UD , 



d^U-tA':U + tUA'^,{A9) 



f;, = z[D^^, Dl] = d^Al - d^Al - i[Al, All 

with r ~ L, R. The pion weak-decay constant and the 
quark condensate in the chiral limit read 



fB„ = -{qq) = 



while the chiral coefficients are 



(47r)^ 



U = 

L5 = 

Lr = 

Ls = 

L9 = 

Hi = 

H2 = 



-2L2 = 
Lg = 0, 



-4Li 



Po 



(47r)2 2Bo' 

_Nc 1_ (_p\_ 

(47r)2 2Nf \2Bo 



(47r)2 
— 2Lio 



A. 

4i?2 



4Bo 
Nc Po 
(4^)2 3 ' 



(47r)2 6 



Po 
12 

_ Po 
24 



Nc P'o 
(47r)2 6 ' 

Nc fp^ 
(4^)2 {Bi 



P'l 
2Bo 



Po 
12 



(50) 
(51) 

(52) 
(53) 
(54) 

(55) 

(56) 

(57) 

(58) 

(59) 



where Np = 2, 3. As we can see, the coefficients 
Li, L2, L^, L^, Lq, Lg, Lio are pure numbers, and coincide 
for convergent integrals with those e xpe cted in the limit 
where the regularization is removed {2l|. The argument 



^ The value of L7 displayed here corresponds to the SU(3) model 
only. For the U(3) model one gets L7 = but then the term of 
Eq. should be added, and the value of L7 is changed. 



anticipating this result in Ref. |4l| has to do with the di- 
mensionless character of the low energy couplings which 
thus involve the zeroth moment po — 1. Note that this 
remarkable result holds without removing the regulariza- 
tion 2 . The fact that Hi is proportional to corresponds 
to a scale-dependent or divergent gauge-field wave func- 
tion, and was observed already in Ref. 'ZT]. Hence, the 
finite piece of Hi depends on the regularization scheme. 

We can use / and L5 in order to determine L7, Lg, Bq 
and H2, which immediately yields 



Nc 



L5 



2Nf SSAn^Nf 



l^R — — 



Nc 



P 



3847r2 64^2 



Ho 



L5 + 



Nc 



P 



1927r2 4^2 



-0.35 • 10-^ 
~ 1.05 • 10-^ 
2.1 • 10"^. 



(60) 



The numerical values displayed here have been obtained 
with the large- TVc value of L5 from Table I. 



B. Curvature contributions 

The curvature contributions to the chiral Lagrangian 
can be written in the form proposed in Ref. and are 
given by 



(61) 



and 



jr{i,R} ^ -LiiRiD^'U^Df.U) - Li2R^''{D^U^D,U) 
- Li3R{x^U + U^x)+H3R'' + HiR^,R^"' 

+ i/5i?M.a/3i?'^''"''. (62) 

Here R^^^, Rfii, and R are the Riemann curvature tensor, 
the Ricci tensor, and the curvature scalar, respectively"^, 

dA o pA o pA 1 pA pa pA pa 

Rtiv = R^tiXv ; R = g^^Rfiv ■ (63) 



The Christoffel symbols are specified in Eq. (jA22p . The 
curvature terms reflect the composite nature of the pseu- 
doscalar fields, since in the considered model they corre- 
spond to the coupling of the gravitational external field 



Actually, the kinetic energy term obtained in Ref. j2]| within 
the ^eta-function regularization was scale dependent, so dimen- 
sional transmutation sets in. If dimensional regularization is 
used, it would lead to a 1/e-divergence, which after renormal- 
ization would also lead to dimensional transmutation. The point 
of the spectral regularization is that dimensional transmutation 
is precluded thanks to the spectral conditions, Eqs. J^, and any 
choice of the spectral function yields the same finite result. 
Note the opposite sign of our definition for the Riemann tensor 
as compared to Ref. We follow Ref. l58l (see Appendix 1X1. 
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at the quark level. After some algebra we get 



where 



Ho = 

Ll3 — 

H3 = 

Hi = 

= 



4 12' 

-2iii = 



Nc Pq 
(47r)2 6 ' 



2 19R„ 6 



(47r)2 12Bo 

(47r)2 ^90' 
(47r)2 720 



(64) 
(65) 

(66) 

(67) 

(68) 

(69) 



Note that there is a finite strong renornialization to New- 
ton's gravitational constant G, since the classical Ein- 
stein's Lagrangian is £ = — i?/(167rG). This correction, 
proportional to the ratio of the hadronic to the Planck 
scale f^Gn/i, is numerically tiny. 



C. Energy-Momentum tensor 

Using the action of Eq. ^ one can compute the energy 
momentum tensor as a functional derivative of the action 
with respect to an external space-time-dependent metric, 
9tiu{x), around the flat space-time metric rj^^i, (we take 
the signature (H ) ), 



1 



r''(x) = 



5T 



(70) 



_i^^d^p(^)(^|{OM.,(iD)-i}|:r), 



where 



Of'' ^ (^MQ-^ + Ydl^) - gt^'^ (i^ _ uj) 



(71) 



In the flat space-time limit, g'^" = 77'^", the chiral La- 
grangean contains only metric contributions and takes 
the form given in Ref. 3, 4], 



£ = +£(2) +£(4) 



(72) 



where 



£(2) ^ £(2,S) 
£(4) ^ £(4,g) 



(73) 



If we do a derivative expansion (see Appendix for de- 
tails) the effective chiral energy-momentum tensor up to 
and including fourth order corrections in the chiral count- 
ing reads ^.5] 



^?(0) + 0(2) + 0(4) 



fiiy 



fiiy 



fiiy 



(74) 



0iO) 
0(2) 
0(4) 



(0) 



(75) 
(76) 



- 2Ln {g^,d^ - d^d,) {D^U^ D^U^) 

- 2L,3{g^.d^~d^d,){x^U + U^x) 

- L12 {g^j^agi^ad^ + gtii^dadp - g^ad^dp - g„adf,dp) 

(77) 



Note that the coefficients Li — Lio appear in 

£(4) 

given by 

Eq. (|48(l . The terms containing Ln — L13 cannot be ob- 
tained by computing the energy momentum tensor from 
the chiral effective Lagrangean in fiat-space time (|72() and 
form this viewpoint are genuine quark contributions to 
6"^" in this model. Actually, the difference between com- 
puting the energy-momentum tensor from an action at 
the quark, i.e. starting from Eq. (|71() . or at the meson 
level, i.e. starting from Eq. (|72l) is 



5T 



5S9 



(78) 



with and denoting the metric and curvature con- 
tributions to the action, is precisely related to the curva- 
ture terms corresponding to the couplings Lii,Li2, and 

^13- 



VI. RESULTS FOR THE MESON DOMINANCE 
MODEL 

The Meson Dominance Model (MDM), developped in 
Ref. [ij, offers a particularly simple realization of the 
SQM and provides an explicit form for the spectral func- 
tion. The quark propagator becomes 



S{p) 



where 



Ps{^) 



pv{uj)i> + Ps{^^)^ 



Zip' 



i> - M(p2 



1 1 1 
2^^(1 -4u;2/M2)5/2' 

1 12P3 
27ri M|(l-4tj2 /Af 1)5/2 ■ 



(79) 

(80) 
(81) 



The vector spectral function, pv(uj)^ is determined by 
imposing vector meson dominance of the pion electro- 
magnetic form factor, from which the identity 



247r2 



(82) 
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is deduced. This relation is subject to chiral correc- 
tions. It is remarkable that such a simple relation pro- 
duces a mass of My = 826MeV for = 93MeV which 
agrees with the value recently obtained in Ref. [s^ . 
With this value of / one gets a vacuum energy of i? = 
-SNpf/Nc - (202-217MeV)'' for Np = 3. In contrast 
to pv{uj), the expression for the scalar spectral function, 
Ps{^), is an educated guess which fulfills the odd spectral 
conditions pi = ps ~ . . . = and reproduces the value 
of the P3 log-moment. The preferred value for the vector 
mass is 



My = 770MeV, 



(83) 



which corresponds to the p-meson mass and which is used 
in the subsequent numerical analysis. 

The integration contour C used in the MDM encircles 
the branch cuts, i.e., starts at —oo + iO, goes around the 
branch point at —Mv/'2, and returns to — oo — iO, with 
the other section obtained by a reflexion with respect to 
the origin ^ij. These two sections are connected with 
semicircles at infinity. The mass function becomes 



M(p2) 
M(0) 



( Ml \ 



5/2 



M2 



( Ml \ 



5/2 



where the constituent quark mass is** 



Mr 



M{0) = - 



5MlNc 



(84) 



(85) 



When Af(p2) ^ p2 ^i-^^^ ^■(p^) = q, such that the 
quark propagator has no poles in the complex p^- 
plane. Instead, it has a cut starting at the branch- 
point = My /A. The exponents reproduce accurately 
the l/(— p^)^/^ behaviour in the deep- Euclidean domain. 
This behavior was seen in the recent QCD lattice simu- 
lation in the Landau gauge, linearly extrapolated to the 
chiral limit jl^. A fit to the data yields 

Mq = 303 ± 24 MeV, 
Ms = 970 ±21 MeV, 

with the optimum value of per degree of freedom equal 
to 0.72, yielding and impressive agreement of M(p^) up 
to = — 16GeV^. Although Z{p^) is not nearly as good 
(c/. Ref. |4^). leaving room for improvement, we think 
it worthwhile to pursue the pattern of chiral symme- 
try breaking which arises in this particular realization 
of the SQM. Incidentally, let us note that if the results 
of Sect, mil are used we get 



47r2 



dp% 



5/2 ' 



(87) 



* In Ref. |47l| there were typographical errors in Eq. (10.6) and 
(10.9), which should carry an extra factor of 2 on the RHS. 



TABLE I: The dimensionless low energy constants (multi- 
plied by 10^) compared with some reference values and other 
models. The errors for SQM in the MDM realization reflect 
the errors in Ms and Mq of Eq. 18611 . 



X ID 


bQiVl 


y^Lir 1 


T f AT b 

Large iVc 


In JL 


Dual 




(iVlUIVl ) 








Large Nc 


L\ 


0.79 


0.53±0.25 


0.9 


0.96 


0.79 




1.58 


0.71±0.27 


1.8 


1.95 


1.58 


L-i 


-3.17 


-2.72±1.12 


-4.3 


-5.21 


-3.17 


Li 

















Ls 


2.0±0.1 


0.91±0.15 


2.1 


1.5 


3.17 


Le 

















Lr 


-0.07±0.01 


-0.32±0.15 


-0.3 






Ls 


0.08±0.04 


0.62±0.20 


0.8 


0.8 


1.18 


L9 


6.33 


5.93±0.43 


7.1 


6.7 


6.33 


Lio 


-3.17 


-4.40±0.70" 


-5.4 


-5.5 


-4.75 


Lii 


1.58 


1.85±0.90^ 


1.6 






Ll2 


-3.17 


-2.7 


-2.7 






Ll3 


0.33± 0.01 


1.7 ±0.80^ 


1.1 







"The two -loop calculation of Ref. |^ 

''Ref. |3 

'^Ref. 13 

''Sec footnote □ 

■^Rcf. 10, ej 

-fRef. 



which reproduces Eq. (|82|l and shows the consistency of 
the approach. For the Meson Dominance model we get 



, MD 
Pi 



I MD 
92 



NcMl 



P3 



■ MD 



Nr. 



5MqMI 

Ml 

6 ' 
5MqMI 

12M2 



Using these values we get 



L5 
Lj 

Ls 



Nr. 



Ml 



967r2 Ml ' 




(89) 
(90) 
(91) 



In the SU(3) case we display our results in Tabled We 
note that the predictions for .^1.2,3,4,6,9,10 £^re common to 
the scheme of ReL ,21^. The values of ^5,7,8 are specific 
both to the SQM and the MDM realization. 

In the SU(2) case we have, with the help of the rela- 
tions given in Ref. [3|, to pass form SU(3) to SU(2) |3]. 
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In the absence of meson loop corrections ^ 



h = 
h = 



-h - --^k - -^k - 



-Nc, 



3 

3M| 



167VcM^o 
75M2 M| ' 



(92) 
(93) 

(94) 



The vector and scalar pion radii are given by Q 



1 



3 



(95) 



respectively. While the vector pion mean squared ra- 
dius reproduces the built-in vector meson dominance 
of the pion e.m. form factor, the scalar radius shows 
that the scalar mass obtained by a fit to the lattice 
quark mass function does correspond to the mass of a 
scalar meson dominating the scalar form factor, {r'^)^'^ = 
0.50±0.01fm. 

The scalar (spin-0) and tensor (spin- 2) component of 
the gravitational form factors, 6*0 and 62 respectively, 
produce the same mean squared radii. 



(96) 



regardless of the particular realization of the spectral 
model. If we saturate the form factors with scalar and 
tensor mesons, /o and /2, we get for their masses 



M 



fo 



Mf^ = 47r/^ VS/TVc = 1105 - 1168MeV. (97) 



depending whether we take / = 88 or 93MeV, respec- 
tively. The experimental value for the lowest tensor me- 
son is Mf^'P- = 1270MeV. As discussed in Ref. ;|, the 
9q (corresponding to the trace of the energy-momentum 
tensor) form factor couples to scalars, whereas the O2 
(corresponding to the traceless combination of 0^^) form 
factor couples to tensor (spin-2) mesons. 

One message is clear from the present model: the 
scalar meson of mass M/g which dominates the energy- 
momentum tensor does not necessarily coincide with the 
scalar meson of mass Ms, which dominates the scalar 
form factor. Actually we have Mfg — \/2Mv, whereas 
Ms is a free quantity. This is natural in the spectral 
approach, where in the chiral limit the scalar form fac- 
tor Fs involves the odd spectral moments, whereas 
involves the even spectral moments. In particular, the 



5 The relations are h = 192-rT^{2Li_ + L3), h = 192tt^L2, h 
2567r2(2L4 +_L5 - ALg - 2Ls), h = 64tt^{2L4 + L5), h 
-1927r2Lio, le = 1927r2L9, hi = 192tt^Lxx , ^13 = 2b&TT'^h 
The constant Z12 is not renormaUzed by the pion loop. 



corresponding mean squared radii are proportional to p'l 
and to poi respectively. Finally, we note the numerical 
value of ^3 = 4.65 obtained in MDM amounts to a shift 
of the pion mass by less than 1% and an increase of 
yielding 89MeV as compared to / = 87MeV. 



VII. THE LARGE-iVc LIMIT AND DUALITY 

Given the fact that our result corresponds to a one- 
quark-loop approximation, we cannot expect our model 
to be better than the leading large-TVc contribution to 
the low-energy parameters, which is made of infinitely 
many resonance exchanges On the other hand, the 
evaluation of these large- iVc contributions requires addi- 
tional, not necessarily unreasonable, assumptions such as 
the convergence of an infinite set of states, and moreover, 
an estimate of the contributions of higher resonances. In 
practice, one works in the Single Resonance Approxima- 
tion (SRA) yielding a reduction of parameters [slfisll. 



2^SRA 



'SRA 



SRA 



1 

= -L 



SRA 



1l 



SRA 



P 



8M2 ' 



(98) 



SRA 



P 



AMV 



L 



SRA 



SRA 



2L 



SRA 



13 



-SRA 
^12 



= -3L 



3LSRA 



SRA 



-SRA 



4M2 ' 

JO 



2A//? 

J2 



(99) 
(100) 
(101) 

(102) 



where /, My and Ms should stand for the leading large- 
Nc contributions to those quantities. To obtain the for- 
mulas for Li till Lio, the pseudoscalar and axial me- 
son contributions have been fine tuned to satisfy the 
VV-AA and SS-PP two point correlation functions high- 
energy-behavior chiral sum rules plus some well converg- 
ing high energy properties of hadronic form factors. (In 
particular, Mp/Ms = Ma/ My = a/2, where Mp is 
the mass of the excited pion. Obviously, more short- 
distance constraints require more resonances. The values 
of Lii,i2,i3 are obtained from the single scalar and ten- 
sor resonance exchange On the one hand, a tensor 
meson is needed in order to provide a non- vanishing L12 
as a minimal hadronic ansatz, on the other hand ten- 
sor mesons do contribute also other LEC's [U, which 
is not taken into account in Eq. (|102(l . Thus, to sim- 
plify the discussion, in what follows we restrict ourselves 
to the non-gravitational couplings Li till Liq. In prac- 
tice, the phenomenological success is achieved by using 
the physical values of the parameters. Note that al- 
though there is predictive power, it is done in terms of 
two dimensionless ratios, f /My and f /Ms- Obviously, 
in the chiral limit we expect both My and Ms to scale 
with /tt. Therefore, in order to preserve the large- A''c 
counting rules one should have My = cyp/y/Nc and 
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= cs/,r/-\/7Vc with cy and eg denoting some Nc- 
independent coefficients. Remarkably, in the SQM the 
low energy parameters depend on two dimensionless ra- 
tios, p'l/ B[) and p'2/B^. It is therefore tempting to deter- 
mine the spectral log-moments from large Nc arguments, 
in a model-independent way. Actually in the Single Res- 
onance Approximation (SRA) we note that the ratios 
Li : L2 : Lg of the SQM agree with those of the SRA. 
The values of L5 and Lg can then be used to determine 
p[, and p'2, respectively, yielding 



/ SRA 

Pi 



P2 



, SRA 



NcMl ' 



6 ' 



(103) 
(104) 



in agreement with Eqs. H89|l and Eq. H82|l . This is not sur- 
prising since the physics of the meson dominance version 
of the SQM and the SRA approximation is alike. The 
only difference is that one cannot deduce from Eqs. H104|) 
the value of the constituent quark mass Mq = Af(0), 
which is given by the ratio of two negative moments 
Mq = /5-i//0-2, Eq- ©■ To determine Mq would re- 
quire computing terms of 0{p^) in the chiral Lagrangian 
and comparing to the SRA at large Nc- 

One can see that it is not possible to match Lg nor 
LiQ. The disagreement with the large- A^c values of 
and iio has to do with the fact that the SS-PP sum 
rule and VV-AA second Weinberg sum rule are violated 
in the present as well as other quark model calculations 
[12, (except for the non-local models, see 0, H^). 
This calls for a modification of our model. The disagree- 
ment has to do with the absence of axial-meson exchange 
in LiQ (1/4 of the total contribution) and pseudoscalar 
meson exchange in Lg (1/4 of the total contribution). 
On the other hand, for the value of / obtained from 
Eq. the constants Li, L2, ^4, ^5, ^6, Ln repro- 
duce the large- A/'c constraints obtained in Ref. [lj|. This 
agreement is confirmed in Table if one corrects for the 
factor 24:n^f^/NcMy = 1.15. One could force L3 to agree 
to the large- A/'c estimate by taking My — Ms- This 
agrees with the observation of the Chiral Unitary ap- 
proach of Ref. [S^; in the large- A^c limit the scalar and 
vector mesons become degenerate ^. Thus, the marriage 
of large- A^c m SRA approximation with our chiral quark 
model calculation produces degenerate scalar and vec- 
tor mesons. Degenerated scalar and vector mesons were 
suggested very early in the context of superconver- 
gent sum rules and have been interpreted more recently 
on the basis of mended symmetries [67l |. Experimental 
claims have been raised [t^ and contested Ffll . 

Direct experimental tests have also been suggested [7^ . 



For Nc = 3,10,20,40 Ref. ^ obtains Ms /My = 
0.58, 0.84, 0.96, 0.98 respectively with Ms and My the real part 
of the poles in the second Riemann sheet. We thank J.R. Pelaez 
for providing these numbers. 



It is clear that whatever sensible modification of the 
SQM is considered, it will only affect Lg and Lio, keep- 
ing the remaining L's. We leave the explicit construction 
of such a modified model for a separate study. Regard- 
less on the particular way to achieve this, we may antici- 
pate already on the consequences for the large Nc in the 
single resonance approximation of taking AIs — My = 
2nf y/6/Nc, yielding the following duality relations. 



2Li 



L2 - 
Nc 



1 



-L, 



1 



-Ls 



3 



rLg 



1 



:Lg 



1 



-Li 



(105) 



1927r2' 

This also implies the set of mass dual relations. 

Ma = Mp = V2Mv = V2Ms = ^TTy/3/NcU ■ (106) 

The new relation M^ = Mp agrees with the experimental 
number within the expected 30% of the large- A'c limit. 
Using Eqs. H95() we obtain 



(r 



2\l/2 



(107) 



These relations are subject to higher 1/A^c and m^r cor- 
rections. We may account for the latter by allowing f-^ to 
vary between the physical value and the value in the chi- 
ral limit. This yields, (r^)^^ = {r^)\/^ = 0.58-0.64 fm. 
The value of the scalar radius is compatible with the one 
obtained in ChPT to two loop P|, 0.78 fm. Going to the 
SU(2) case, in the dual large- A'c model we get 

- h ^l2 = \h = \h = \h = \h = Nc, (108) 

whereas the recently extracted values obtained at the two 
loop level from the analysis of tttt scattering 1^ and vector 
and scalar form factors at the two loop level are 

= -0.4 ± 0.6 , I2 = 6.0 ± 1.3 , I3 = 2.9 ± 2.4, 
[4 = 4.4 ± 0.2 , [5 = 13.0 ± 1.0, fe = 16.0 ± 1.0. 

(109) 

The I coefficients are in a sense more suitable for compar- 
ison with ChPT since the chiral loop generates a constant 
shift in all of them by the same amount, c = log(/i^/m^). 
Thus, it makes sense to compare the differences where 
chiral logs are canceled. We find 



h 
I3 

k 



2Nc 



I = ^ 
^ 3 

- ^ 5Nc 

^ 3 

h = 4.Nc 

h = 5Nc 



(Exp. 6.4 ± 1.4), 
(Exp. 3.3 ±2.4), 

(Exp. 4.8 ±0.4), 

(Exp. 13.4 ± 1.1), 
(Exp. 16.4 ± 1.1), 



(110) 



where the errors have been added in quadrature. As 
we can see, the agreement is excellent, within the un- 
certainties, and suggests accuracy of the order of 1/N^ 
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rather than the standard a priori l/7Vc error estimate. 
The constant pion loop shift can be accommodated with 
a scale ^ = 513 ± 200MeV, comparable to the p meson 
mass. Taking Eqs. I|1U2|I . corresponding to the SRA with 
the physical values / = 93MeV_^/s = lOOOMeV, and 
Mv = 770MeY, as_done in Ref. [ll, yields_Z2 - = 8.3, 
h-h= 6.2, h-h= 6.2, ^5 - ^1 = 15.2, h-li = 18.7. 
More reasonable values are obtained by taking Ms — 
eOOMeV, but then the SRA relation Mp = V^Ms pre- 
dicts a too low value of the excited pion state. The 
present discussion favours phenomenologically the dual 
relations (|105|l as compared to the SRA relations H1U2|) 
with physical parameters. 



VIII. CONCLUSIONS 

In the present work we have studied the chiral ex- 
pansion of the recently proposed Spectral Quark Model 
in the presence of electroweak and gravitational exter- 
nal sources. The model is based on a Lehman repre- 
sentation for the quark propagator with an unconven- 
tional spectral function, which is genuinely a complex 
function with cuts in terms of the spectral mass. We 
have written down the effective action which reproduces 
the Ward-Takahashi identities presented in the previous 
work. Thanks to an infinite set of spectral conditions 
demanded from the power like factorization property of 
form factors at high energies, we have been able to show 
that the corresponding chiral anomalous contribution to 
the action is properly normalized without removing the 
regularization. Moreover, the non-anomalous contribu- 
tion to the action can be written in the long wavelength 
limit in terms of 13 low energy constants. The numerical 
values are in reasonable agreement with the phenomeno- 
logical expectations, although some discrepancies do oc- 
cur for Ls and Liq. In some cases they can be naturally 
explained as failures in reproducing some chiral short dis- 
tance constraints which suggest that the model needs to 
be improved. On the other hand, if one tries to match 
the remaining non-gravitational LEC's to large Nc pre- 
dictions in the single resonance approximation, a further 
reduction of parameters takes place. In particular, one 
finds the best agreement for degenerate scalar and vector 
mesons. 

We have estimated for the first time in the frame- 
work of chiral quark models the gravitational LEC's Ln, 
Li2 and ii3, describing the coupling to external grav- 
itational sources. These LEC's depend on curvature 
properties of the curved space-time metric. This calcu- 
lation allows a determination of some matrix elements 
of the energy momentum tensor. Our analysis suggests 
that the scalar meson coupling to the quark conden- 
sate moqq and the scalar meson coupling to the trace 
of the energy momentum tensor 6*^ do not necessarily 
coincide. Clearly, these two operators behave differently 
under chiral symmetry, since moqq vanishes in the chi- 
ral limit whereas 0^ does not. This point is in itself 



rather intriguing and deserves further investigation. We 
note here that this fact materializes in our model be- 
cause these two scalar mesons depend on odd and even 
spectral moments respectively. On the other hand, we 
obtain M/„ = M/^ = V^My = V^Ms = a 
very reasonable result if we take into account the large 
Nc nature of the one quark loop approximation. Fur- 
ther quark-meson duality relations have been discussed, 
allowing a rather successful! determination of the best 
known LEC's, consistent up to the experimental errors 
with the best known values up to two loop accuracy. 
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APPENDIX A: DERIVATIVE EXPANSION AND 
USEFUL IDENTITIES 

1. Reduction to a vector like theory and 
transformation properties 

The Dirac operator can be rewritten as 

D = -DrPr + -DlPl, (A1) 
with the projection operators on parity 

^fl. = ^(l + 75), Pl = \{1-15), (A2) 
such that for a Dirac spinor one has 

^fl = PR'f, *L = Pi*. (A3) 
The right and left Dirac operators are given by 

iDl = i^ + 4^-M\ (A4) 

with 

M = s + ip + ujU , = s - ip + ujU\ {A5) 

= t;/- + ^ Al^v^"- a". (A6) 

the quark mass matrix is included in the scalar field s. 
Under left-right unitary transformations, and ^Ir, one 
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has the following properties, 

■^R -> n^-^R, *L -> (A7) 

u ^lU^^r, c/t ^^J^c/t^7t_, (as) 
^ nRA^i.ii'j^ + inRd^'n^R, (A9) 
^ + (aio) 

The chiral covariant derivatives and field strength tensors 



F. 



fiiy 



r^R,L. (All) 
behave as follows under local chiral transformations: 



D^-^R ^rD^^r, 
D^^L ^ ^lD^-^r, 



(A12) 
(A13) 
(A14) 
(A15) 



2. Coupling the Spectral Quark Model to Gravity 

The coupling of fcrmions to gravity is well known (see, 
6.17., Ref. ,73]) but not in the context of chiral quark 
models. We review it here for completeness and to fix 
our notation. We use the tetrad formalism of curved 
space-time (for conventions see, e. 5., Ref. 58]). Given the 
metric tensor we get a local basis of orthogonal vectors 
(tetrads or vierbein). 



g'^-(a;)=e^(.x)e^(x)ry 



AB 



(A16) 



with rj^^ — diag(l, —1, —1, —1) for a flat Minkowski met- 
ric. These vectors fulfill the orthogonality relations. 



6it 



M A 



t-A uu A A fJ. 



(A17) 



Under the coordinate — > x'^{x) and frame x^ — s- 
AgX^ transformations the transformation properties of 
the tetrad are 



frame transformations according to Eq. ljA18|l as follows, 
The covariant derivative is defined as 

"'tJ.^vA — "t^^vA^ ^ v^i^\A + ^ iiX^vA+^ABfiJ^u ' 

(A21) 

where the Riemann connection is given by the Christoffel 
symbols, 

- ^5"" {5a5m- + 5m.9a. - d.g^x] , (A22) 

which are symmetric in the lower indices, F^^ = rj^_^ 
(we assume here no torsion). In order to preserve the 
covariance of the tetrad mapping we must have 

dp^e^,A = d^_,eA,u - T^^e^.A + t^ABf^e^ = 0. (A23) 

In addition, the condition d^g^'' = 0, implying 

df^rjAB (^ABfj^ + (^BAfi = 0, (A24) 

requires an antisymmetric spin connection, ujABf^ — 
-^BAi,, given by 



(A25) 



The frame and coordinate covariant derivative is de- 
fined according to the spin of the corresponding field. For 
a spin-0 C/, spin-1/2, ^E", spin-1, A^, and spin 3/2, ^ 
fields the transformation properties are 

U{x) ^ U{x), 

^{x) S{K{x))^{x), (A26) 
A,{x) ^A^x), (A27) 

- |j5(A(x))*.(x). (A28) 

For infinitesimal Lorentz transformations A;^ = '^b + 
with eAB — ~^BA one has 5(A) = 1 — jOABf^^^ with 
(JAB defined below (see Eq. ljA34|) '). 

For a scalar (spin-0) field we have the standard defini- 
tion 



dx'" 
dxt^ 



e-; ^ A^(x)e^, 



(A18) 



respectively. The tetrads map coordinate tensors into 
frame tensors (which transform covariantly under local 
Lorentz transformations), for instance 



(A19) 



Frame tensors are invariant under coordinate transfor- 
mations x'^. For a general tensor T"^ greek in- 
dices transform covariantly under coordinate transfor- 
mations while latin indices transform covariantly under 



For a vector (spin-1), one has 



dpiAi, — dfj^Ajj — T^^Ax, 



fulfilling the property 



[d^, d^] Aa ~ R\ij.uAx 



(A29) 
(A30) 
(A31) 



with the Riemann curvature tensor given by Eq. H63|l . 
The coordinate and Lorentz covariant derivative for 
Dirac fermions (spin 1/2) is defined as 



(A32) 
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where lo^ is the Cartan spin connection, 



-CT UJABfi, 



and 



(A33) 



(A34) 



with the "fA are fixed a;— independent Dirac matrices (we 
use the conventions of Ref. [7j|) fulfiUing the standard 
flat space anticommutation rules, 

7^7^ + 7^7-^ = 277-4^. (A35) 

The space-time dependent Dirac matrices are defined as 

j^{x)^-fAe^{x) (A36) 

and fulfill 

-/^'{x)-f•'ix) + 7''(a;)7^(a;) = 2g^''(x). (A37) 

The covariant derivative of a frame (x-independent) 
Dirac matrix (behaving as the adjoint representation 
^•^) is 

^^7^ = dp.jA - i [t^/i, 1a] + tOABfj^l^ = 0. (A38) 

Thus, we obtain the useful identity for the coordinate 
(and cc-dependent) Dirac matrix, 



d^^-fy{x) = 0, 



(A39) 



which implies that for the free Dirac operator the order is 
irrelevant = ^^^{x)d^'^ = d^'y^(x)'^ . For a mixed ten- 
sor (spin-3/2) the frame and coordinate covariant deriva- 
tive reads 



(A40) 



Applying the previous definition to d^^'^ one gets the use- 
ful formulas 



[d^,,dy] * 



fT"''i?a/3M.*, (A41) 

= {(a^ - lu^) [V^g^'" (a. - loj,)] vf} , 

(A42) 



a „/3 „AB 



is an antisymmetric x-dependent 



where cr"'^ 
matrix. 

Gauge fields can be included by the standard minimal 
substitution rule yielding the covariant derivative for a 
fermion. 



V^* = [d^ - iA^) ^. 



(A43) 



With this notation the full Dirac operator in the presence 
of external vector, axial-vector, scalar, pseudoscalar and 
gravitational fields reads as in Eq. lO, where 



4 = 7^(a;H,(x), 



(A44) 



and the pseudoscalar Dirac matrix in the curved case is 
defined as 



75(2;) 



1 



4!, 



£^''"''7^(^)7. (^)7a(x)7/3(x) 



^e^^'^^7A7B7c7D = 75- 



(A45) 



Here g{x) = Aet{g^'') since det(e^)2 = det(g'^'') with 
g0i23 _ 2 (both in the frame as well as in the coordinate 
sense) . 

The full coordinate, frame and chiral gauge covariant 
derivative for pseudoscalar (spin-0), Dirac spinor (spin- 
1/2) and a Rarita-Schwinger spinor (spin 3/2) fields are 
given by the following formulas, 



V^U ^ d^U - i[v^,U]- i{a^,U}, 



(A46) 

and they correspond to replacing the derivative by the 
frame and coordinate covariant derivative, d^, in 

the chiral covariant derivative Z?^. Note that with this 
definition neither V^V^'^ ^ V^V^^E' nor Vf^V^U are co- 
ordinate covariant since the second derivative does not 
include the Riemann connection F^,,. 



3. The Second order Operator 

In the absence of gravitational sources, the normal par- 
ity contribution can be obtained from the second order 
operator (see Eq. 

D5D = [pl + iM^pL-ipnM^ +M^M]Pr 

(A47) 

Gravitational fields can be coupled by covariantizing first 
the Dirac operator, i.e. making 9^ — s- or — s- 
and taking into account that since a spinor field is a co- 
ordinate scalar we have 



(A48) 



The same reasoning can be applied to the coordinate 
scalar y^*, yielding 



(A49) 



This means that we can assume Pl.r = ^ l.r. when act- 
ing on spinor field as follows 

D5D* = [yl+iMfL~'yfRM+M'^M]PR:^ 

+ [y'i + iM^fL-iyRM'^+MM'^]PL-^. 

(A50) 
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If we include the gauge fields we have two vector like 
theories with left and right gauge fields and re- 
spectively. Suppressing momentarily the left and right 
labels we have 



(A51) 



where the use of the identity 

= [D^, D,] ^ + '-ct^PR^p^,-^ (A52) 

has been made. The coordinate and frame invariant 
Laplacian for a Dirac spinor is given by 



1 



-g 



(A53) 



Note that for a Dirac spinor field 5* the operator I?^ 
contains the spin connection. Reinserting the right and 
left chiral notation the second order operator takes the 
suitable form 



D.D 



1 



[V,{V^g^'-'V,)]+V, (A54) 



with 



any derivative 9^ first order. This implies in particular 
that R^'^"^, R^'^ and R are taken to be of second order. 
Finally, the external scalar and pseudoscalar fields s and 
p are taken to be second order as well. Thus, the mul- 
tiplicative operator V — is at least first order in the 
chiral counting. To the computed order 0{p'^) in the heat 
kernel expansion one has to go up to 04. The contribu- 
tions can be separated into the flat space non-vanishing 
contributions and the curvature contributions generated 
by quantum effects. Using the form suggested in [t^I we 
have 



0-2 



= 1, 



cj'-V+-R, 
6 



^— R n R*"""/' 1— 7? 7?^'' -I- — f T 

180 1 Qo /^"^ ^ 1" 



1 9 1 9 
6 



180 
1 



as = 
where 



30 
1 
6 
1 

24 



12 



:R 



6 



(A59) 



and 



Vl 



V = VrPr + VlPl 



(A55) 



(A56) 



4. Derivative expansion 



MM 



v^v = v^v„v. 



(A60) 
(A6I) 



Clearly, the heat kernel coefficients depend on the spec- 
tral mass a; in a polynomial fashion. Using the integrals 



— (it) 



^ — 2g — ircj 



[ioyriz - 2) (A62) 



we get for integer z = n and after using the spectral 
conditions, Eq. Q, the normal parity contribution of the 
action takes the form 



We use the proper-time representation, 

Tr log (D5D) = -Tr / —e'"^'^ + C, (A57) 
^0 ^ 

with C and infinite constant. The form of the opera- 
tor D5D in Eq. (jA54p is suitable to make a heat kernel 
expansion in curved space-time as the one of Ref. [75j |. 
For a review see e.g. 76] and references therein. In our 
particular case, before undertaking the heat kernel ex- 
pansion we separate a uP' contribution from the operator 
D5D which we treat exactly. 



x) (A58) 



(47rir)2 



n=0 



a2n{x) (ir)" . 



The derivative expansion is done by considering U zeroth 
order the vector and axial fields and first order, and 



d^x^/^ i dujp{uj) 



-ATrlogDsD = 

2 ^ ^ 2(47r2) 

X tr(— icj'' logw^ao -f cj^ logw^fli 
- log(tj^/^^)a2 + + \ai + . . .) 

d^x^g (£(2) + £(4) + 

(A63) 

After evaluation of the Dirac traces, the second order 
Lagrangean is 



£(2) = 



(47r) 



p(cj)| -cj2 logtj2(v^[/tvA'[/) 



+ 2uj^ XogLu^ {m^U + U'^m) + log ^ (i?) | , 

(A64) 
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whereas the fourth order becomes 

Nr. 



£(4) 



(47r)2 
1 
6 



+ i:\ogu;'{{Fp' + {F^^r) 



1 o 1 
TD^ I 

144 90" 



+ ^((V^C/V,C/t)2)-i((V,[/V^C/t)2) 

+ logw^w2^2(mtm) + ((mtC/ + C/V)2)) 

- iw(V^C/^V^?7(m^C/ + C/^m)) 

- logw2w(V^[/^V'^m + V^mty^f/) 

- ujloguj^^R{U''m + m^'U) + ^RVf,U^V^'U)y 

(A65) 

Note that up to this order the moments pa ~ 1, pi — 
and p2 = as well as the log-moments Pq, p\ and p'2 
appear. 



5. Equations of Motion 



We define, 



X = 2Bom = 2Ba {s + ip) . 



(A66) 



For on-shell pseudoscalars one may minimize the action 
at lowest order 



5(2) ^ Z! 



d*3 



X (V^[/tV^[/+(xt[/ + C/tx) 
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R), 



(A67) 



to obtain the equations of motion. Since U is unitary, 
U^U = 1, we have that the variations on U and C/^ are 
not independent of each other, SU^U + WSU = 0. For 
SU(3)-flavour one has, in addition, to impose the con- 
dition DetJ7 = 1. One can treat U and C/^ indepen- 
dently by introducing a term in the Lagrangian of the 
form {KWU — iAlog/7) where the Lagrange multipliers 
are A, a hermitean matrix, and A, a real c-number. Thus, 
the EOM are 

V^C/ = x+(A-iA)C/, 
V^t/I' = + U\k + i\), (A68) 



where 



V^C/ = -^D^ (V^g'^-'D^U) . (A69) 



Combining these two equations, we get 

[/tV^?7- V^C/tf/ = U^x-x''U-2iX. (A70) 

Taking the trace and using the condition that for a ma- 
trix with DetU = 1 one has {U^W ^U) = and hence 
{WV^U - V^WU) = 0, we get 



6i 



(A71) 



thus 



On the other hand A is given by 

2A = W^U^U + UW^U^ - (xC/^ + x^t^)- 



(A72) 



(A73) 



Using the identities deduced form the unitarity condition 

Wu = 1, 

t/^V^C/ + V^J/tf/ = (A74) 
C/^V2f7 + V^C/tf/ = ~2Vf,U^V''U, (A75) 



and combining them with the previous Eqs. (jA73IA75p 
we get the identities 



+ l(xtc/-c/tx)' 



(A76) 



and 



(xty^cz + v^c/tx) = 2{x^x)-l{{x^u + u^xf) 
- ((xtc/ + c/tx)v^c/tv^c/) 

+ lix^U + U^xf- (A77) 
6 

In the case of the U{3) group one has DetU = e^"^"^^ ^ 
1 and the last two terms involving (x^C/ ± C/^x)^ in 
Eqs. HA76I) and (|A77|) should be dropped. (See the dis- 
cussion before Eq. H59|) ) The result can be further sim- 
plified using the integral identity 



-g{\J^'\J''U'^Vi'\/''U) ^ j d^Xy/^{\7^U^\7^U) 
f^R^^iWU^W'U), (A78) 



which can be deduced from Eq. (|A3ip applied to V^t/. 
Finally, we also have the SU(3) identity 



((V^l/^V,[/)2) = -2(V^C/tv^C/)' 



(A79) 



+ (V^i7tv,;7)2 + -(V^C/tv^[/)2. 
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Once the identities l|T76|) . l|T77|l . (|X78|l and l|J79|) have 
been used one can make the substitute the coordinate- 
frame-covariant derivative by the covariant derivative, 
i.e., V^U — D'^U, since the pseudoscalar matrix U is 
a coordinate and frame scalar. In that way Eqs. (|48|l and 
(|()^ are deduced. 

In four dimensions, one can reduce the form of the 
curvature contributions to the Lagrangean if the Gauss- 
Bonnet theorem is used in Eq. H62|l . namely that 

is a topological invariant (the Eulcr number ) and hence 

= (A81) 

under metric deformations, — > -I- 5g^^. This rela- 
tion was not taken into account in Ref. Q but it does not 
affect the calculation of the energy momentum tensor in 
flat space, Eq. (f7H|l . 



Derivative expansion for first order differential 
operators 



tor. The derivative expansion of the Dirac operator can 
be done using the identity 



'ip - M~ujU' 



1 



(27r)4^ + i^-7W-wC/' 

(A82) 



where the differential operator acts on the right. This 
formula can be justified by requiring vector gauge invari- 
ance of the action or by using the asymmetric ver- 
sion of the Wigner transformation presented in Ref. |55| | . 
Expanding in powers of D and Ai and squaring the de- 
nominator we get 



1 



d-^k 



fc2 



n+1 



^7 {2tiY 

•1=0 ^ ' 

ujU^) [{ip - M) + loU^)]" . 

(A83) 



As we see the definition of the action involves the Dirac 
operator D only, which is a first order differential opera- 



In this way Eq. 1441) can be derived. 
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